


New Results in the Calculation of Modulation Products 

By W. R. BENNETT 

A new method of computing modulation products by means of multiple 
Fourier series is described. The method is used to obtain for the problem of 
modulation of a two-frequency wave by a rectifier a solution which is con- 
siderably simpler than any hitherto known. 

THE problem of computing modulation products has long been 
recognized as being of fundamental importance in communication 
engineering. Heretofore certain quite fundamental modulation prob- 
lems have been attacked by methods which are difficult to justify from 
the standpoint of mathematical rigor and some of the solutions ob- 
tained have been in the form of complicated infinite series that are not 
easy to use in practical computations. In this paper these problems 
are solved by means of a new method which is mathematically sound 
and whicli yields results in a form well suited for purposes of com- 
putation. 

The analysis here given applies specifically to the case of two fre- 
quencies applied to a modulator of the "cut off " type ; i.e., a modulator 
which operates by virtue of its being insensitive to input changes 
throughout a particular range of values. A simple rectifying charac- 
teristic forms a convenient basis of approximation for study of such 
modulators, and hence we consider in detail methods of calculating 
modulation in rectifiers when two frequencies are applied. Applica- 
tions to certain other types of modulation problems and to the case of 
more than two applied frequencies are discussed briefly at the close. 

Half Wave Linear Rectifier — Two Applied Frequencies 
We shall define a half wave linear rectifier as a device which delivers 
no output when the applied voltage is negative and delivers an output 
wave proportional to the applied voltage when the applied voltage is 
positive. We may take the constant of proportionality as unity since 
its only effect is to multiply the entire solution by a constant. Assume 
the input voltage e(t) to be specified by 

e(t) = P cos (pt + d p ) + Q cos (qt + $,). (1) 

The output wave will then consist of the positive lobes of the above 
function with the negative lobes replaced by zero intervals. It is 
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convenient to represent the amplitude ratio Q/P by k, and without 
loss of generality to take 

P > and < k ^ 1. (2) 

The problem we now consider is the resolution of the output wave into 
sinusoidal waves, a complete solution requiring the determination of 
the frequencies present, their amplitudes, and their phase relations. 

The method of solution used employs the auxiliary function of two 
independent variables f(x, y) denned by 

f(x, y) = P(cos x 4- k cos y), cos x + k cos y S: 0, 1 
= 0, cos x -f- k cos y < 0. J 

It is clear that the function /(.v, y) may be represented by a surface 
which does not pass below the ry-plane and which coincides with the 
.ry-plane throughout certain regions which are bounded by the multi- 
branched curve, , 

cos x -f k cos y = 0. (4) 

If either x or y is increased or decreased by any multiple of 2w, the 
value of f(x, y) is unchanged. Hence /(.r, 3/) is a periodic function of 
x and y, and if its value is known for every point in the rectangle 
bounded by y = ± t, x = ± t say, the value of the function may be 
determined for any point in the entire .ry-plane. 

From the above considerations we are led to investigate the expan- 
sion of f(x, y) in a double Fourier series in x and y. We may readily 
verify that the function satisfies any one of several sets of sufficient 
conditions ' to make such an expansion valid. We may write the 
expansion thus: 

f(x, y) = £ £ \_A± mn cos (mx ± ny) + B± mn sin (nix ± ny)\ (5) 

m=0 n=0 

with the summation to be extended over both the upper and lower of 
the ambiguous signs except when m or n is zero, in which case one value 
only is taken (it is immaterial which one) ; when m and n are both zero, 
we divide the coefficient Am by two in order that all the i4 -coefficients 
may be expressed by the same formula. Determining the coefficients 
by the usual method of multiplying both sides of (5) by the factor the 
coefficient of which is to be found and integrating both sides throughout 
the rectangle bounded by x = ±71-, y = ±tt, we obtain: 
1 Hobson, "Theory of Functions of a Real Variable," Vol. 2, p. 710. 
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A± m „ = -5-5- I I fix, y) cos (m* ± ny)dydx, 

B± mn =-^~t I I /(*» y) sin ( w * ± n y) d y dx - 



(6) 



We now return to our original problem of representing the positive 
lobes of a two-frequency wave as a sum of sinusoidal components. We 
may apply the double Fourier series expansion of f(x, y), which must 
hold for all values of x and y, to the special case in which x and y are 
linear functions of the time. If we let 

x = pt + dpA ,js 

y = qt + 9 ,\ 

the function f(x, y) represents the rectified two-frequency wave as a 
function of time. The values of x and y which are used lie on the 
straight line, 

y=^ + e g -p pt (8) 

which is obtained by eliminating / from (7). A representation of 
f(x, y) valid for the entire xy-plane must of course hold for values of 
x and y on this .straight line. Hence we may substitute the values of 
x and y given by (7) directly into the double Fourier series (5), and 
the result will evidently be an expression for the rectifier output in 
terms of discrete frequencies of the type (mp ± nq)/2ir. The phase 
angle of the typical component is md p ± nd q and the amplitude is 
expressed by (6). 

The solution is thereby reduced to the evaluation of the definite 
double integrals of (6). Three different methods of reducing these 
integrals have been investigated, and it appears that each has certain 
peculiar advantages and points of interest. We shall consider them 
separately. 

I. Straightforward Geometric Method 

In this method, which yields remarkably simple results in a direct 
manner, we determine the boundaries of the region throughout which 
f(x, y) vanishes and substitute appropriate limits in the integrals to 
exclude this region from the area of integration. When this exclusion 
has been accomplished, f(x, y) may be replaced in the integral by 
cos x + k cos y. The boundary between zero and non-zero values of 
f(x, y) is the curve (4), which has two branches crossing the rectangle 
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over which the integration is performed. The non-zero values of 
f(x, y) lie in the shaded region of Fig. 1. From the symmetry of the 
region about the x and y axes we deduce at once that the sine coeffi- 
cients, B± mn , must vanish and that the cosine coefficients, A± mn , may 
be obtained by integrating throughout one quadrant only and multi- 



y o 




Fig. 1 — Region of integration for the determination of the coefficients in the double 
Fourier series expansion of/ (x, y). 



plying by four. We therefore obtain, on substitution of the proper 
limits, 

IP r 
in = — T I cos nyay 

1 Jo 



A mn = A, 



•arc cob ( — *• oos y} 

> I (cos x + k cos y) cos mxdx (9) 



Jo 
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This expression gives the amplitude of the typical component of 
frequency {mp ± nq) /2v. The remaining steps are concerned merely 
with the calculation of the integral (9) for particular values of m and n. 

It will suffice to work through one example in detail and give the 
results in tabular form for the other products up to the fourth order. 
The second order side frequencies, (p ± q)/2ir, will be taken as a 
typical case. 

By direct substitution 

9 T3 /»t /»arc coa ( — A" cos y) 

An - -j ) cos ydy I (cos x + k cos y) cos xdx. (10) 

Performing the inner integration and substituting the limits for x, we 
obtain : 

A n = — I cos y[arc cos {—k cosy) + fccosy Vl — k 2 cos 2 y~]dy. (11) 
* Jo 

Considering separately the integral, 

cos y arc cos (— k cos 7)^7, 






integrate once by parts, letting 



u = arc cos (— k cos y), 
dv = cos ydy. 



The result is, after combining with the remainder of the integral for .An, 

. kP C v sin 2 y + cos 2 y(l — k 2 cos 2 y) , , r9 s 

7T J V 1 — & 2 cos'' y 



Now substituting 

cos y = z, 
we obtain 



tt 2 J >/(i-a«)(l- W) 

This is a standard elliptic form. 2 It is convenient here to let 

2 It may be remarked that a large number of the integrals required in the evaluation 
of the coefficients are listed by D. Biererts de Haan, "Nouvelles Tables d'Integrales 
Definies." See in particular Tables 8 and 12, pages 34 and 39. 
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(I z m 
dz. (14) 
V(l - z 2 )(l - kh 2 ) ■ 



By differentiating the expression z m_3 V(l — z 2 )(l — k 2 z 2 ), we may 
easily derive the useful recurrence formula: 

7 _ (m - 2)(1 + F)Z m _ 2 - (m - 3)Z m - 4 M _. 

Zm " (m - \)k 2 (15) 

We may now calculate the value of Z m for even values of m in terms of 
Z and Z 2 . Z is a complete elliptic integral of the first kind which we 
shall designate as usual by K; i.e., 

Z = K = f dz = H \ 1 - k 2 sin 2 0dO. (16) 

J \(1-S 2 )(1 - kW) X 



Furthermore from. the identity: 

z 2 = i r i _ /i - kh 

\ (1 - S 2 )(l - PP) & L\(l -S 2 )(l - «) \ 1 - 2 2 

we have 

1 



, (17) 



Z 2 =^(K-E), (18) 

where £ is a complete elliptic integral of the second kind defined by 

E= I J 1 ~ k *f dz = \ 2 \ 1 - k 2 sin 2 Odd. (19) 

Now making use of (15), we calculate Z 4 in terms of Z 2 and Z and 
get finally: 

„ (2 + £ 2 ^ - 2(1 + £ 2 )E 

Zj - Jkl ( 20 ) 

We can then evaluate (13) in terms of K and E. The result is 

An - ^ [(1 + * 2 )£ - (1 - fc 2 )*]. (21) 

The process of evaluating the other coefficients is quite similar. 
Results are listed in Table I. 

Convenient tables of K and E may be found in Peirce's Short Table 
of Integrals (page 121), Byerly's Integral Calculus, and the Jahnke 
und Emde tables. For a very extensive set of tables, see Legendre's 
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Traits des Fonctions Elliptiques. Numerical calculation of the coef- 
ficients making use of these tables and the formulae listed in Table I 
is a quite simple process. Curves of the coefficients as functions of k 
have been calculated in this way and are plotted in Fig. 2. 




'. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 Q 0.8 

RATIO OF FUNDAMENTAL AMPLITUDES , -p- = K 

Fig. 2 — Curves showing amplitudes of modulation products in output of half wave 
linear rectifier when input wave consists of two frequencies. 

It is perhaps worth noting that the special case of equal funda- 
mental amplitudes (P = Q or k = 1) yields the simple result, 






A = 



8(-)"P 

\_{m + w) 2 - l][(w - w) 2 - 1> 2 ' 



(22) 



where m -f- n is even. When m + n is odd and greater than one, 
A mn is zero. 
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II. Fourier Series Method 

The second method is of interest because it obtains the same results 
as the only previously known solution, 3 which is in terms of infinite 
series involving Bessel functions. The fact that the results agree is a 
check on the validity of certain doubtful rearrangements of multiple 
series necessary in the process by which these results were originally 
obtained. Furthermore by comparison with the corresponding results 
of the first method we can sum the infinite series in terms of complete 
elliptic integrals; a number of interesting mathematical theorems are 
thus proved, which have been made the basis of a paper by the author 
in the December, 1932 issue of the Bulletin of the American Mathe- 
matical Society. 

By expanding the function: 



tt 

4>(u) = — j> — c < w < 



-» < u < c 



(23) 



in a Fourier series in u, we may verify that: 




= u, < u < c 



If we let u = P cos x + Q cos y, the left hand member of (24) is equal 
to f(x, y) provided \P\ + \Q\ < c. With this restriction on c, we 
may substitute the resulting expression for f(x, y) in the integrand of 
(6), and no change in the limits of integration are required. Term by 
term integration of the series may be justified without difficulty, and 
making use of well known definite integrals, we obtain finally: 

- m+n+2 _ J m ( TrP ) Jn ( " vQ ) 

a -**(-) • y __V i / V £ L, (25) 

mn ir^ ' it! (2r-.l) a 

where m + n is an even integer. When m + n = 0, the extra term 
e/4 must be added. When m + n is odd and greater than one, the 
value of A mn is zero; when m + n = 1, the values are Am = P/2, 
A n = Q/2. 

Peterson and Keith obtained the above result 3 by substituting 

3 Peterson and Keith, "Grid Current Modulation," Belt System Technical Journal, 
Vol. 7, pp. 138-9, January, 1928. 
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u = P cos x + Q cos y in the left hand member of (24), applying 
Jacobi's expansions in series of Bessel coefficients, and rearranging the 
resulting triple series. It appears that it is much more difficult to 
justify the series rearrangement than term by term integration. From 
the results obtained by the first method it follows that the series in 
(25), which might be termed a generalized Schlomilch series, 4 is sum- 
mable in terms of elliptic integrals. 

III. Trigonometric Integral Method 

Following a suggestion of Mr. S. O. Rice, we may make use of the 
following relation: 



u u C sin u\ j. 

o + - I — x — d\ = u, u > 

2 7T J X 

= 0, u < O.J 



(26) 



Evidently if we substitute u = P cos x + Q cos y, the left hand member 
of (26) represents the function f(x, y) and may be substituted in the 
integrand of (6) without change in the limits. Interchange of the 
order of integration may then be justified without difficulty and the 
following result is obtained in terms of a special case of the integral of 
Weber and Schafheitlin: 

where m + n is even and greater than zero. When m + n = 0, the 
above integral should be replaced by an infinite contour integral taken 
along the real axis except for an indentation to avoid the origin and 
with all other quantities remaining the same except for a division by 
two. For all even order modulation products it may now be deduced ■ 
that: 



(_)— + ' r f "» + "-! ]k , P 



A m „ = 



2,r(« + i)r( '"-;' + 3 ) 



The case of m + n = requires a special investigation, which shows 
that (28) holds for this case also. 

4 Cf. Watson, "Theory of Bessel Functions," Chapter XIX. 
6 Watson, "Theory of Bessel Functions," p. 401. 
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The hypergeometric function in (28) may always be expressed in 
terms of K and E by successive applications of recurrence formulae and 
use of the known relations: 



■K 



K-inr v u 



IT 



1 1 



E-iF[-v *; 1; 



2 1 \ 2 r 



(29) 



By means of the hypergeometric recurrence formulae we may also 
show that 

. 2[(m — l)k 2 + n— \~\Am-\, „-i + (m + n — 5)fei4 m -a, n-2 ,- m 

when w + n is even. A discussion of the hypergeometric function and 
a derivation of (30) are given in the appendix. 

From (30), we can compute successively all even order modulation 
products starting with say Aoo and An known. If negative subscripts 
occur in applying the formula, they may be replaced by positive sub- 
scripts without changing the validity of the results; this is proved in 
the appendix. 

Half Wave Square Law Rectifier — Two Applied Frequencies 
The solution for two frequencies applied to a square law rectifier, or 
in fact to any rectifier operating on an integer power law, can be 
obtained in a manner quite similar to that used in solving the linear 
rectifier. In the case of a square law rectifier, we have to represent 
the function 

f(x, y) = P 2 (cos x + k cos y) 2 , cos x + k cos y > 1 . . 
= 0, cos x + k cos y < 0. J 

Going through the same steps with this function that we did with that 
of (3), we find that the amplitudes of the modulation products can be 
expressed in terms of K and E as in the case of the linear rectifier; the 
results are listed in Table I. A set of curves is plotted in Fig. 3. 

We may also show that 

m+n+l c , „ Jm ( " irP\j n [ irQj 

A = (-) 2 — T ^ - '- - - - (32) 

A mn - { ) ^ i, (2f _ 1)3 IJz; 

when m + n is odd and greater than one and c > \P\ + \Q\- For 
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Fig. 3 — Curves showing amplitudes of modulation products in output of half wave 
square law rectifier when input wave consists of two frequencies. 
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Av> and Am. we must add cP/2 and cQ/2 respectively. The value of 
A mn is zero for m + n even and greater than two; the other even order 
products are listed in Table I. Another form of the result for odd order 
products is 

7T Jo A 



or 



k «p> r fm + n-2 



«m» — \ 



2*r(« + i)r ( 



m — n + 4 
~2 



xF ^n±n- lf n-m-2 . % + 1; ^ (3f) 
A three term recurrence formula for odd order products is: 

2[(W— l)fe + W— l]i4 w -l. w -l+(ffl+W— 6)M w - 8 . n-2 /..x 

4 »- " (m+n+2)* ' ^ ; 

When P = Q, and w + w is odd, 

A - 64 (-) -+ip» } 

Amn ~ (m 2 - K 2 )[(m + w) 2 - 4][(m - «) 2 - 4> 2 ^ ; 

Other Applications and Results 
The solution for any full wave rectifier can be obtained from the 
solution for the corresponding half wave rectifier. Thus we may easily 
show that the output of a full wave linear rectifier contains neither of 
the fundamentals and that the amplitudes of all other modulation 
products are twice as large as the corresponding amplitudes in the 
output of a half wave linear rectifier. It is also evident that by super- 
posing the solutions for the linear and square law rectifiers we can 
obtain the solution for a quadratic law rectifier having an output equal 
to aie(t) + a 2 [e(/)] 2 when e(t) is positive and no output when e{t) is 
negative. Biased rectifiers, peak choppers, and saturating devices 
can be solved by the same methods used above, the solution of course 
becoming more complicated for the more complicated kinds of char- 
acteristics. Nor is the method restricted to " cut off " type modula- 
tion. Curvature type modulators can be treated in the same way and 
in many cases solution by the above method is simpler than by the 
usual power series expansion. The method also appears to have 
promise in the solution of magnetic modulation problems, where the 
effect of hysteresis must be considered. 
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When three frequencies are applied, a triple Fourier series is required, 
and in the general case of n frequencies, a Fourier series in n variables 
would be used. The work becomes more complicated as the number of 
frequencies increases, but there is no theoretical limitation. 

In conclusion the writer wishes to express his appreciation of the 
valuable advice of Messrs. T. C. Fry and L. A. MacColl on the tech- 
nical features of the paper. 

Appendix 

The hypergeometric function F(a, /3; 7; 2) may be denned by the 
power series : 

lly 217(7 + 1) 

When any one of the three quantities a, /3, 7 is increased or decreased 
by unity a new hypergeometric function is formed which is said to be 
contiguous to the first. Gauss listed fifteen linear relations which 
connect F(a, /3; 7; 2) with pairs of its contiguous functions. In deriv- 
ing the recurrence formula for A mn we require difference relations 
between functions which are not contiguous, but the required relations 
may be obtained from those listed by Gauss by a process of substitution 
and elimination. 

We shall find it convenient to designate F(a, /3; 7; 2) by F, F(a + 1, 
/3; 7; 2) by F a+ , F(a + 1, /3; 7 — 1; z) by F a+y _, etc., and to let 

m + n — 3 n— m — \ 
at = 2 ; P = 2 ' 7 = n ' Z = 

In this notation, Equation (27) becomes: 



\—2~+ 1 ^ (m -\- n — \ 



m+n 

The corresponding expressions for A m -\, n -i and ^4,„_2, »-a are by direct 
substitution: 



■n-t-n , , _> 

(-)~r ( '" + » ~ 3 ) k-'P 



A m -i, „_i — ■ : — rr .F, 

2xr( w )r /, "- w - j ^ 



m+n 



(-)~"'r ( w+ 2 "- 5 ) k-*P 

A m— 2, n— 2 — — J j ^T r a —y—. 

27iT(» - i)r I J 
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Thus a recurrence relation expressing A mn in terms of A m -i,n-i and 
A m -2, „-2 evidently requires a relation between F a+y+ , F, and F a - y — 
Referring to Gauss' tables, 6 we find 

(7 - a - i)F + aF a+ - (7 - 1)^- = 0, 
7(1 - z)F- yF a . + (7 - P)zF a+ = 0. 

From the second of these two equations we form two more equations by 
substituting a + 1 for a in one case and 7 — 1 for 7 in the other, 
giving 

7 (1 - z)F a+ - yF+ (7 - P)zF a+y+ = 

(7 - 1)(1 - s)F 7 _ - (7 - 1)^-7- + (7 " 1 - P)zF - 0. 

Now eliminating F a+ and F y - from the first, third, and fourth of the 
equations, we obtain 

a(|8 - 7)s-Fa +7+ + 7 [T - 1 + (« - fi)*]F 
- 7(7 - l)F a . y - = 0, 

which is the relation desired. Substituting the value of F a+y+ in 
terms of A mn , F in terms of A m -i, „_i, and F a _ 7 _ in terms of A m - 2l „_ 2 
gives the recurrence formula of Equation (30). 

In using (30) we may find, as for instance in calculating A m0 , A m \, 
A 0n , Am, that the right hand member involves coefficients with nega- 
tive subscripts. A simple rule for treating such cases may be demon- 
strated as follows. We first note that if we replace m by — m in 
(28) the value of the right hand member is unchanged. 7 Hence since 
(30) is derivable directly from (28), we conclude that correct results 
are obtained from (30) if we adopt the convention, 

Si— m, n **-mn- 

The case of n negative is a little more difficult because if n is a 

negative integer in (28), an indeterminate form results. However, 

making use of the result just obtained on the interchangeability of sign 

of the subscripts, m, m — 1, m — 2 in (30), we can demonstrate a 

8 Gauss, Werke, Bd. Ill, page 130. The equations used here are numbered (5) 
and (8) by Gauss. 

t If we express (-)-/»r ( W + \ - ' ) /v ( " = \ ± 3 ) in terms of (-)-« 

Y |~ ffl + n ~ 1 j I y (^^- — n — J by successive applications of the recurrence 

formula for the gamma function, we find the two quantities are equivalent. Chang- 
ing the sign of m in the hypergeometric function merely interchanges a and /3, and 
hence does not change the value of the function. 
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similar rule for the subscripts n,n— 1, n — 2, valid when (30) is used. 
For example, by direct application of (30), we deduce that 

A _ 2[(1 - m)k 2 + n + l]^i- m , n+1 + (n - m - l)kA. m . „ 

A 2 ~ m ' n+2 " (n - m + 5)* 

Now since it is known that we may replace the subscripts 2 — m, 
1 — m, and —m by m — 2, m — 1, and m respectively, we may show 
that 

a _ 2[(/» — \)k 2 — n — l]i4 w _i, B +i + (w — w + 5)fe^ m _ 2 , n+2 

(w - n + 1)* 

which is exactly equivalent to the relation we get if we replace n by — n 
throughout in (30) and then substitute A m , _„ = A mn , A m , _„_i = 

Am, n+]» n«i, — n— 2 = Am, n+2- 

It may be remarked that it would be incorrect to base a proof of 
interchangeability of sign of subscripts on (27) because the equivalence 
of (27) and (28) has not been demonstrated for a sufficient range of 
values of m and n. 



